arXiv:1509.00050v2 [hep-ph] 13 Mar 2016 


New Prospects in Fixed Target Searches for Dark Forces 
with the SeaQuest Experiment at Fermilab 

S. Gardner* 

Department of Physies and Astronomy, University of Kentucky, Lexington, Kentucky 40506-0055 USA 

R. J. Holtl^ 

Physics Division, Argonne National Laboratory, Argonne, IL 60439 USA 

A. S. Tadepalli^ 

Department of Physics and Astronomy, Rutgers University, Piseataway, NJ 08854 USA 

An intense, 120 GeV proton beam incident on an extremely long, iron target generates enormous 
numbers of light-mass particles that also decay within that target. If one of these particles decays 
to a final state with a hidden gauge boson, or if such a particle Is produced as a result of the 
initial collision, then that weakly interacting, hidden-sector particle may traverse the remainder of 
the target and be detected downstream through its possible decay to an e''"e“, fi~, or 
hnal state. These conditions can be realized through an extension of the SeaQuest experiment at 
Fermilab, and in this initial investigation we consider how it can serve as an ultrasensitive probe of 
hidden vector gauge forces, both Abelian and non-Abelian. A light, weakly coupled hidden sector 
may well explain the dark matter established through astrophysical observations, and the proposed 
search can provide tangible evidence for its existence — or, alternatively, constrain a “sea” of 
possibilities. 
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I. INTRODUCTION 

Searches for new physics have long been motivated 
by the seeming inadequacies of the Standard Model 
(SM). Some are theoretical and motivate searches for 
new physics at high-energy colliders, such as those that 
would help explain the origin of the weak scale, v = 
« 174GeV [1-3]. Others include its in¬ 
abilities to explain either dark matter or dark energy 
and their relative predominance over visible matter in 
the cosmic energy budget, as deduced from astrometric 
observations [4-6]. The missing new physics can appear 
either at high energies, at short distance scales [7-10], 
or at low energies, at long distance scales [11-18]. Solu¬ 
tions to the dark matter problem could conceivably come 
from either source [19-23]. New, long-distance effects are 
both possible and discoverable if the new, light degrees 
of freedom couple to SM fields in a weak but yet appre¬ 
ciable way. Such operators, or portals, to possible hidden 
sectors have been discussed extensively in the context of 
apparent cosmic and gamma ray anomalies, see Refs. [24- 
27], e.g., and have been proposed as an explanation of 
the muon g — 2 anomaly [28]. (We also note earlier work 
in which a light U(l) gauge boson directly couples to 
SM fermions to explain the muon g — 2 [29-31], as well 
as earlier astrophysical [32], anomalies.) In this paper 
we discuss the discovery prospects of an ultrasensitive, 
broad-band search for new, long-distance physics, made 
possible through an extension of the SeaQuest E906 ex¬ 
periment at Fermilab. 

The success of the SM in describing known particle 
phenomenology motivates a framework in which new 


physics appears as additions to the SM through effective 
operators Oi of mass dimension four or higher. The as¬ 
sociated coupling constants are characterized by CijtG 
where n > 0, Ci is dimensionless, and A is the energy 
scale of new physics. As we have noted, these additions 
are thought to be either at high energy scales for which 
A > V [33-35] with G ~ 0(1) or at low energy scales for 
which A ^ n or n = 0 with G ^ 0(1). In the latter cat¬ 
egory, the most effort has been invested in operators for 
which n = 0, in part because their appearance does not 
usually require the inclusion of additional new physics 
to be theoretically consistent at high energy scales [23]. 
Dark photon searches fall into this class. Higher-mass- 
dimension portals are also possible, but have received 
much less attention. Their coupling to SM particles is 
expected to be much smaller since n > 0. The experi¬ 
ment we consider, in which very small couplings can still 
be appreciable, thus serves as an ideal hunting ground 
for such effects. 

Much experimental effort has been invested in searches 
for rare exotic particles through so-called beam-dump ex¬ 
periments, in which detectors are mounted downstream 
of a particle beam stopped in a target. Our current in¬ 
vestigation concerns the discovery prospects of an exper¬ 
iment of this class, and a schematic is shown in Fig. 1. As 
emphasized by Ref. [37], the potential parameter space 
for a dark photon — or, indeed, for any particle probed 
via a n = 0 connector — is vast, in both candidate mass 
and mixing angle. Beam-dump experiments that involve 
displaced vertices for dark photon production and decay 
are largely sensitive to small mixing angles, with electron 
and proton beam-dump experiments giving comparable 
constraints [37]. Dedicated efforts have been made to 
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FIG. 1. Schematic of the SeaQuest spectrometer layout [36]. 
The 120 GeV proton beam from the Fermilab Main Injec¬ 
tor approaches the spectrometer through a 25-cm-long hole 
of 2.5 cm in diameter in the 5-m-long solid iron magnet. An 
A' generated in the first meter of the beam dump traverses 
the Focusing Magnet (FMAG) without being affected by the 
magnetic field and can decay in the fiducial region into a lep¬ 
ton pair, or a pion pair (upon upgrade). Stations 1, 2, and 3 
comprise a series of drift chambers and an array of scintillator 
hodoscope paddles used for track reconstruction and trigger¬ 
ing purposes. The 3-m-long air-gap KTeV Magnet (KMAG) 
is used to focus the muons back into the spectrometer to facil¬ 
itate momentum measurements. The 1-m long iron absorber 
wall is followed by an array of proportional tubes used for 
muon identification. 



M^, (MeV) 


address the remaining holes in parameter space [37-39], 
with many recently completed searches and reanalyses of 
earlier ones [37, 40-55] with many more proposed and 
under development [56-65]. 

The extension of the SeaQuest experiment we consider 
can contribute to this effort in different ways. Not only 
can it probe new regions of dark-photon parameter space, 
leading either to a dark-photon discovery or a refinement 
of that phase space, but it can also be used to probe dark 
forces that enter solely through their mixing with QCD 
degrees of freedom. In this latter case, proton beam- 
dump experiments play a special role, particularly if the 
downstream spectrometer can detect pions. In what fol¬ 
lows we expound on the discovery prospects of the experi¬ 
ment shown in Fig. 1. For reference, the manner in which 
the existing SeaQuest spectrometer can contribute to a 
dark photon search, with reference to efforts under devel¬ 
opment worldwide, is shown in Fig. 2. We note that the 
SeaQuest experiment can probe part of the dark-photon 
parameter space to be probed by the SHiP experiment at 
CERN [64]. We also consider the discovery prospects as¬ 
sociated with a spectrometer upgrade to permit electron 
and pion detection. 

We now sketch the content of the sections to follow. 
We begin with an overview of hidden portal models, high¬ 
lighting, in particular, the various ways in which quark 
and gluon degrees of freedom can also play a role. We 
then proceed to describe the specihc manner in which 
hidden portals can be probed at SeaQuest. In this ini¬ 
tial investigation we place a particular focus on the ra¬ 
diative decays of the light mesons 7r° and t]. We ex¬ 
pect light mesons to be produced copiously in a proton 
beam dump experiment [75], and their radiative decays 


FIG. 2. Plot shows the projection contours of the coupling 
constant e as a function of dark photon mass for four dif¬ 
ferent processes that could be used to search for dark photons 
at SeaQuest. Regions I and II are bounded by the contour 
plots for rj —> yA' —>■ 7 e''"e“ and rj —> yA' —>■ , respec¬ 

tively, whereas regions III and IV refer to the limits inferred 
from use of the proton bremsstrahlung production mecha¬ 
nism, followed by A' —>■ e^e“ (III) and A' —>■ (IV) de¬ 

cay. The area excluded by electron beam dump experiments 
E137 ]37, 40], E141 [41], and the searches by BABAR ]50], 
GHARM [43, 44], NA48/2 [47] are bounded by solid lines at 90 
% GL, whereas those excluded by u-Cal I (7r°)[51] and u-Cal 
I (p-Brem)[52] are bounded by solid lines at 95 % CL. Also, 
the planned sensitivities of APEX (full run) [56], HPS [57], 
DarkLight [58] (all at 90% CL) and LHCb]65] (at 95% CL) 
are shown as dotted lines for comparison. We omit the an¬ 
ticipated limits from VEPP-3 [59], Refs. [61, 62], Mu3e ]63], 
and MESA [60], which all probe lighter masses, as well as 
Ref. [64], for visual clarity. The region above e = 10~® (not 
shown in the figure) has been excluded by several experiments 
such as E774 [42], APEX (test run) [45], HADES [66], KLOE 
[46], PHENIX [48], MAMI [49], along with the 2cr exclusion 
limit obtained from {g — 2) ^ [67]. Approximate limits at still 
weaker mixing angles from the LSND experiment ]38, 68, 69] 
and from astrophysical considerations ]37, 70-74] have been 
omitted. Note that the limits shown all assume that decays 
of the A' to the invisible sector are nonexistent. 


are controlled by the chiral anomaly even if the final 
state contains strongly interacting particles. The proton 
bremsstrahlung contour for SeaQuest in Fig. 2 was pro¬ 
duced by following the method outlined in Ref. [52] and 
the simulation techniques described in sec. IV.A. Finally, 
we turn to a discussion of the experimental prospects, 
illustrating concretely how SeaQuest can probe dark 
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forces, before offering our concluding summary. 

II. HIDDEN SECTOR PORTALS 

The known interactions and particle content of the SM 
are richly complicated, and this itself suggests that the 
bulk of the matter content in the Universe should be 
similarly complex. Existing constraints on its content, 
however, are minimal. Nevertheless, the observation that 
dark matter is stable over at least Gyr time scales begs 
an explanation, and it is natural to think that a gauge 
symmetry of the hidden sector can provide it. The mat¬ 
ter content of such a hidden sector need only interact 
gravitationally with the matter we know — no other in¬ 
teractions are required, though it has become popular 
in recent years to build theoretical models of cosmoge- 
nesis that tie the generation of dark matter with that 
of the cosmic baryon asymmetry, solving two problems 
at once [76-78]. Although the existence of hidden sector 
gauge interactions could potentially impact the morphol¬ 
ogy of dwarf galaxies [79] and have other observational 
consequences [80-83], the ability to discover such hidden 
sectors may rest on the manner they can connect to the 
particles and interactions of the SM. 

It has been popular to consider portals that consist of 
operators that do not require new high-energy physics 
for theoretical consistency, so that n = 0, or less. An 
economical summary of such portals [23, 38] is 

- H^H{AS + XS^) - YnLHN , (1) 

where is the held-strength tensor associated with 
U(1 )y in the SM, H is the SM Higgs field, and A is a 
SM left-handed (lepton) doublet. The explicit hidden 
sector degrees of freedom are a Held-strength tensor 
associated with a hidden U(l) symmetry, a scalar S', and 
a fermion N. These new degrees of freedom can couple to 
further hidden sector degrees of freedom, and the latter 
can be richly diverse. 

The new gauge boson degrees of freedom can be dark 
photons or Z’s, or both, depending on the manner in 
which their mass is generated [67, 84]. If the dark pho¬ 
ton is massless, kinetic mixing with the visible photon 
would engender a DM electric charge [85], and the result¬ 
ing constraints can be severe [79, 86, 87], though they can 
also be weakened by either introducing multiple species 
with a net dark matter electric charge of zero [88] or by 
making the dark-matter particle a (nearly) electrically 
neutral composite [80, 81, 89, 90]. The latter, “dark 
atom” scenarios possess an additional length scale whose 
existence is also constrained by cosmic microwave back¬ 
ground (CMB) observations [91]. In this paper we focus 
on hidden vector gauge bosons that range from some 10 
MeV to 700 MeV in mass, for which the noted astro- 
physical constraints do not operate. We note in pass¬ 
ing, however, that extremely weakly coupled dark pho¬ 
tons with such mass scales are nevertheless constrained 
by supernova cooling [37], as well as CMB and big-bang 


nucleosynthesis (BBN) [72], considerations. A minimal 
enlargement of the SM in the presence of hidden gauge 
bosons is of form [37-39], 

£ = /IsM + V -, (2) 

where V^i, = dfj,A'^ — dvA'^. The kinetic mixing term 
B^''Vf_ii, engenders, upon diagonalization and field redef¬ 
inition, a photon with a A! admixture controlled by the 
small parameter k. Specifically, A'^ —)■ A^ — eA'^, where 
s = KCOs9y\/- In constrast, mixing with the Z is sup¬ 
pressed by a nominally large ratio of masses, namely 
by 0{em\,/M^) [37, 84, 92]. However, with an en¬ 
larged Higgs sector, such as in the two Higgs doublet 
model, mixing with both the photon and Z can occur, 
and the candidate dark gauge boson becomes a dark Z 
or “Zd” [67, 84, 93, 94]. In this case, after diagonalization 
of the kinetic mixing term, the photon and Z acquire a 
small admixture of Zd which couples to SM electromag¬ 
netic and weak-neutral currents as per [67, 84, 93, 94] 

Tldark Z = ~ (^6 ^ J^c)Zdfi, (3) 

Z cos Uw 

with the earlier dark photon model emerging in the Sz = 
0 limit. The appearance oi Z — Zd mixing gives rise to 
low-energy parity violating effects as well [67, 84, 94]. 
Nevertheless, it is apparent that dark photon searches 
also restrict the Zd- 

Thus far we have considered kinetic-mixing models as¬ 
sociated with an Abelian gauge symmetry. The gauge bo¬ 
son mass, be it that of A or can be arranged through 
a Stueckelberg mechanism [95, 96], but it can also be gen¬ 
erated though a SM Higgs mechanism in the hidden sec¬ 
tor [24, 97] ^. The discovery of a massive Abelian gauge 
boson can thus implicitly hint to the existence of non- 
Abelian hidden-sector interactions as well. Non-Abelian 
interactions can also appear explicitly, and we now delve 
into this possibility directly and consider, in particular, 
how QCD degrees of freedom can serve as a portal to a 
possible hidden sector. 

It has long been thought that new matter with QCD- 
like interactions could exist [98-100], though the internal 
color symmetry of QCD would seem to preclude vector 
portals of the sort we have considered thus far. Never¬ 
theless, portals with QCD degrees of freedom need not 
yield higher mass-dimension operators, nor do they nec¬ 
essarily require new UV degrees of freedom for theoreti¬ 
cal consistency. We highlight some of the possibilities in 
what follows. We loosely organize our discussion in terms 
of portals of increasing mass dimension, beginning with 
portals arising from a gauged baryon vector current and 
then turning to portals employing multi-quark operators 
and gluons. 


^ The hidden scalar can also function as a portal to the hidden 
sector, as per Eq. (1). 
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The possibility of a light U(l) gauge boson B asso¬ 
ciated with a gauged baryon vector current is a notion 
of some standing [101-110], though much of its focus has 
concerned the consequences of its kinetic mixing with the 
U(l)v sector of the SM [105-107]. Recently Tulin [109] 
has noted the possibility of a mass-dimension-four con¬ 
nector in terms of quarks, namely, 

>Cb = , (4) 

where is the new gauge held that couples to baryon 
number. Since U(1 )b is anomalous in the SM, its mixing 
with SM helds requires new UV degrees of freedom for 
theoretical consistency [23, 101, 109, 110], though such 
constraints have been successfully implemented, note, 
e.g., Refs. [102-104, 108, 110]. Since the gauge cou¬ 
pling (/B is universal for all quarks, the B is taken to 
be strictly isoscalar, so that B —)• tt+tt” does not occur 
if G parity is not violated [109]. For B bosons in the 
mass window of M^r Mb 620 MeV the primary de¬ 
cay mode is thus B —>■ Tr^y. Mixing of the B with the 
photon can appear through radiative corrections [105- 
107], so that B —> can occur as well. The B 

can be probed through the radiative decays of the light 
mesons [101, 109]. A distinct signature of this particular 
model is the appearance of R —>■ Tr^y decay, a final state 
which the SeaQuest spectrometer cannot detect. Limits 
on the A', however, translate to those on R in a model- 
dependent way [109]. It is worth noting, however, that 
the R need not be strictly isoscalar. As established from 
phenomenological studies of the nucleon-nucleon force, 
both the p and w mesons couple to the baryon vector 
current because the nucleon-nucleon force is not charge 
independent, though the coupling of the lo is roughly a 
factor of 5 larger than that of the p [111]. Generalizing 
Eq. (4) to the form 

Cw = gB'U'^^uB'^ + g%,d-i^dB'^ + ... , (5) 

we term this gauge boson B' and relegate the contribu¬ 
tions of other quarks to the ellipsis. Thus B' —)■ tt+tt” 
can occur without breaking G-parity, as had been as¬ 
sumed in Ref. [101]. The R' can be probed through 
the decay rj —> R'y —)■ 7r+7r“y, which is accessible at 
SeaQuest. Other models can be probed in this man¬ 
ner as well. For example, in the model of Dobrescu and 
Frugiuele [110], the gauge boson that couples to baryon 
number is a leptophobic Z', where the interactions of the 
Z' to quarks, for the first generation, is of form 

>CqZ' = ^ {zqQlI^Ql + ZuURj^UR + ZddRYdR) Z'^ . 

( 6 ) 

The field Qr is a left-handed quark doublet and ur and 
dR are right-handed singlets — we see that the Z' can 
couple to isovector combinations of quarks as well. 

We now turn to the possibility of portals comprised of 
operators of mass dimension greater than four. The exis¬ 
tence of such operators can require the existence of new 


UV physics for theoretical consistency, but this does not 
negate them as a possibility. In this case it is possible to 
connect to non-Abelian hidden sector degrees of freedom 
directly. Starting in dimension six we note the following 
possibilities namely, 

k' 

-j^{gLiR)it.qLiR)){q'L(R)i^q'LiR)) ( 7 ) 

where q' is a hidden sector quark, admitting the possi¬ 
bility of a hidden sector “pion” [38], e.g., for which dark- 
matter models have recently been developed [113, 114]. 
We note Refs. [115, 116] as alternate models of non- 
Abelian hidden sectors. Explicitly non-Abelian portals 
also appear, such as [97] 

f^tr($“R;jtr(|.'R^''‘^), (8) 


where is the QCD Held strength tensor and 

denotes an analogous non-Abelian object, though it need 
not reside in a theory of SU(3) color. Note that the ap¬ 
pearance of new heavy degrees of freedom is also explicit 
through that of the heavy scalars and $ that trans¬ 
form under the adjoint representations of their respective 
groups [97]. Finally, in dimension eight it is also possible 
to have a pure glue-glue connector 

tr(R“'^'^R;jtr(Ri'^R'>'Q, (9) 


which has also served as the basis for “hidden valley” 
models [117] that can possess striking collider signa¬ 
tures [117-120]. This connector is also noteworthy as 
an UV “complete” extension of the SM. That is, oper¬ 
ators with mass dimension greater than four that break 
no SM symmetries do not require additional UV physics 
to be viable, an idea exploited in the dark-matter model 
of Ref. [116]. 

The various higher mass dimension connectors we have 
discussed appear in connection with QCD couplings, so 
that “infrared slavery” at low energy can partially off¬ 
set the ratio of mass scale suppression associated with 
their higher mass dimension. Indeed, at low energies we 
should recast the connectors we have considered in terms 
of hadronic degrees of freedom. A pertinent model is pro¬ 
vided by the hadronic kinetic mixing model of Ref. [121], 
which is based on the hidden local symmetry model of 
QCD [122, 123], in which the p mesons function as ef¬ 
fective gauge bosons of the strong interaction, and on 
vector-meson dominance [124, 125]. That is [121], 


11 E 

r . — /a^i/ I ^ a J a 

'^mix — ^ 

9 2 

7Tl 7Tl / 

+ , (10) 


^ We note, however, that axion degrees of freedom can connect 
to SM fermions in dimension five [38], though their couplings to 
quarks are constrained by BBN [112]. 
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where ^ denotes the hadronic vector current and “ 
are the gauge bosons of a hidden local SU(2) symme¬ 
try with = d^pv^°' — di,p^il^°' [125]. This model 
resembles the dark photon models discussed earlier but 
contains two massive vector fields; its possible footprints 
at low energies, in beta decay, have been considered by 
Gardner and He [121], The kinetic mixing terms can 
be removed through the field redefinition —)• 

thus yielding a coupling of the hadronic vector current to 
p'“. Here the can be probed at SeaQuest. Indeed all 
the quark-level models we have considered that permit 
a hidden-sector contribution to a 7r“'"7r“ final state can 
feed the low-energy constant £ in Eq. (10). In this sense 
a p' search is rather generic, as it constrains the B', lep- 
tophobic Z', and higher-dimension QCD connectors we 
have considered as well. 


III. PROBING HIDDEN SECTORS AT 
SEAQUEST 

Enormous numbers of light mesons would be produced 
in the collision of an intense, 120 GeV proton beam with a 
thick iron target, and we propose to search for the light, 
hidden sector particles that may be produced through 
their decays. Moreover, dark gauge bosons can be pro¬ 
duced directly through initial-state radiation from a high- 
energy proton beam, followed by a hard proton-iron col¬ 
lision; here we employ the formalism of Ref. [52]. 

The models discussed in the previous section can all 
be probed through an extension of SeaQuest, though the 
distinctive features of a dark Z — its ability to medi¬ 
ate parity violation — would not be apparent. Never¬ 
theless, sensitivity to a parity-violating dark gauge bo¬ 
son (Zd) should result from the use of a polarized pro¬ 
ton beam [126] via the initial-state radiation production 
mechanism Specifically, the Z^ can be produced via 
either a vector or axial-vector coupling. These processes 
can interfere to give different i'^£~ yields as a function 
of the beam helicity. 

Irrespective of such developments, searches for hidden 
dark forces in atomic parity violation, as well as in parity- 
violating electron scattering, play an important comple¬ 
mentary role [94]. 

The following decay chains are examples of processes 
that can be studied at SeaQuest. We note that A' 
searches automatically limit the Zd as well. 

• 7r° —)■ yA' —>■ ye+e”: to search for A'. 

• rj —>■ yA' — )■ y£+f“ with £ G e,p,: to search for A'. 

The ability to detect different lepton final states 


® Parity violation could be studied through use of either a polarized 
beam or an exceptionally long, polarized target, but the latter 
is not planned for the extension of the SeaQuest experiment we 
consider here. 


offers the possibility of testing the kinetic mixing 
mechanism. 

• Tj —>■ yA' —>■ y7r+7r“: to search for A' — or p', a 
“strongly” interacting hidden sector gauge boson. 
The presence of the p' would be signalled if no ac¬ 
companying events were observed in £'^£~. 

• proton bremstrahlung of A' with A' -g£+,£- 
and £ e e,/i or A' —)■ tt+tt”: to search for A'. 
Study of A! —)■ tt+tt” also permits a search for a 
leptophobic dark gauge boson Z' [110]. 

• Drell-Yan production of A' with A! ->£+,£“ 
and £ G e,^: to search for A! — studying the 
possibility of A' —> tt+tt” also opens sensitivity to 
a leptophobic dark gauge boson. This perturbative 
QCD mechanism for dark gauge boson production 
is under study [127]. 

Although we focus on the possibility of new vector 
gauge forces in this paper, we wish to emphasize that 
the discovery prospects of the type of search we discuss 
span a much broader horizon. We note that the study 
of light meson decays in this context include these addi¬ 
tional possibilities. 

• rj —>■ TT^hd —> 7r°7r+7r“: to search for a dark Higgs 
hd, noting Refs. [128, 129], as well as for scalar dark 
matter that mixes with QCD gluons [116]. 

• -G TT^Zd ^ with £ G e,/r: to search 

for Zd [94]. 

We now evaluate each of the radiative decay processes 
in turn. 

A. The decay —>■ yA' —>■ ■y£~^£~ 

The decay amplitude 7r° —)• yy(*) is controlled by the 
primitive axial anomaly [130, 131], even if QCD correc¬ 
tions that would alter its form only vanish for on-mass- 
shell photons [132]. Consequently, the partial width for 
Dalitz decay is usually normalized relative to that for 
7r° —> yy [133-135]. To validate our procedures we begin 
by computing the partial width for Dalitz decay at tree 
level in the SM [133-135], which yields, after a partial 
integration over the three-body phase space, 

r(7r°ye+e") = / dx lf(l,0, x)l^ ——— 

61T J^2 X 

x/i(i + d). ,11) 

where Tq = r(7r° yy), x = s/M^o, r = 2me/M^o, 
X, and s is the invariant mass of the e“'"e 
pair. We note that /(l,0,x) is the 7r° —)• yy* transi¬ 
tion form factor (noting /(p^/M^o, fci/Af^o, for 

7r°(p) —> y*(fci)y(fc 2 ) decay [134]), where /(1,0,0) = 1 
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appears in 7 r° —)■ 77 . The transition form factor has 
been studied extensively because of its connection to 
the analysis of the hadronic-light-by-light contribution 
to the anomalous magnetic moment of the muon, g — 2, 
and a model-independent representation, based on Fade 
approximates, of the transition form factor data from 
CELLO, CLEO, BaBar, and Belle, yields 

/(l,0,a;) = 1-I--I--I-e>(x^), (12) 


with = 0.0324(12)stat(19)sys and 6,^ = 

(1.06(9)stat(25)sp) X 10“^ [136], though a disper¬ 

sive approach yields similar results with smaller errors: 

= 0.0307(6) and = (1.10(2)) x lO-^ [137]. For 
reference an analysis of tt^ —)■ 76 “'" e“ in two-flavor chiral 
perturbation theory (ChPT) enlarged by electromag¬ 
netism yields a,r = 0.029(5) [135], whereas ChPT at one 
loop with /X = TOp yields a^n- = 0.036 [138]. Computing 
r( 7 r° —>■ 7 e“''e“)/ro from Eq. (11) using empirical inputs 
from Ref. [139] yields 1.185 x 10“^ if the transition form 
factor is set to one, in agreement with Refs. [133, 135]. 
Including /(l,0,x) using the central values of Ref. [136] 
yields 1.188 x 10 “^, whereas the experimental result 
r( 7 r 0 ^ e+e- 7 )/ro = (1.188(30)) x lO'^ [139]. We 
conclude that our framework works very well. The 
three-body partial width via an A' intermediate state 
follows from inserting the factor 


(^(x - XA'f +XA'Ta') 


(13) 


under the phase space integral, where xa' = rn\,jM^ 
and = Ta'/M-j^, to yield 


pi 

r(7r^ —)■ 7^^ —)■ 7e“^e“) = To / dx —e^|/(l, 0, x)| 

J j -2 7T 

2(1 — x)^x 


2 _ + ■ 


(^(X - XA'f +X^/f^,^ 


(14) 


Existing searches for hidden gauge bosons, such as A! 
and Zdi have presumed that there are no lighter par¬ 
ticles within the hidden sector, so that these particles 
decay only to SM particles. In this case this implies 
B{A'{Zd) —>■ e+e“) = I. Consequently, since Ta' ~ 
ae^ 1 , we can employ the narrow width approxima¬ 
tion, replacing 


(^{X - XA'f + XA'f A') 
to yield 

r(7r° —)■ '-)A! —)■ 7e“''e“) = 


-^ ^_' b (x - XA') ■ 

(15) 

|/(l,0,XA')7(7r° ^ 7^') 
x, 8 (xl'^ e+e"), (16) 


where 


r( 7 r° ^ 7 X 1 ') = 2 e 2 



(17) 


and 


r(xi' 


s^e ) = -«£ mA' 


a-4 


Ml 


1-f 2 



(18) 

We note Pa' = r(xl' —>• e+e ) by assumption. As a 
result, Eq. (16) can be viewed as the product of two se¬ 
quential two-body decays, in which the A! appears on its 
mass shell, moderated by the ir^ transition form factor. 
If, rather, xa' > 1, then the integral in Eq. (14) is regular 
even if Pa' = 0. We set this possibility aside, however, 
both here and in what follows because an on-mass-shell 
particle is needed to evaluate its transit through the iron 
beam stop. It is worth emphasizing that the transition 
form factor we have noted in Eq. (12) is developed for 
the low-momentum transfer regime x ^ 1 , whereas xa' 
can be of 0(1). Using a Fade approximate of P|(s) 
form [136]^, which satisfies the asymptotic limit from 
perturbative QCD [141] for the transition form factor, 
we have checked that for our application its numerical 
effects are so slight that we can neglect it with impunity. 


B. The decay g —>■ 7 A' —>■ 

To evaluate rj —>■ yl'^l~ decay in the SM we need only 
replace and Wg —>■ mi in Eq. (11) and modify 

the transition form factor. A new parametrization [142] 
of the 77 —)■ 7*7 form factor, including the latest mea¬ 
surement of rj —)■ e“''e “7 [143], and employing rational 
approximates [144], has recently become available. With 
X —>■ and as x —>■ 0 (with x > 0) we have 

/^(l, 0, x) = 1 -I- bnX + c^x^ -I- dnX^ + 0{xf , (19) 

where 67 ^ = 0.576(ll)stat(4)sys5 c?? — 0-339(15)stat(5)sys5 
and drj = 0.200(14)stat(18)sys [142]. Since the inte¬ 
gral of Eq. (11) includes x = 1 as well, we also eval¬ 
uate /,,(l, 0 ,x) using the complete rational approxi¬ 
mate^, noting that time-like data is available up to 
s ~ 0.22 GeV^ [142]. Employing Eq. (11) and Eq. (19) 
to compute r (77 —)• 7 e+e“)/r 2 yields 1.67 x 10“^ (and 
1.62 X 10“^ without the form factor) to compare with 
the experimental ratio 1.75(06) x 10“^ [139]. Comput¬ 
ing, rather, T{r] ^ 77 i+/x“)/r(j yields 8.17 x 10“^ (and 
5.51 X 10“^ without the form factor) to compare with 
the experimental ratio 7.87(13) x 10“^ [139]. Using 
the complete rational approximate makes for little im¬ 
pact, changing only the value of r (?7 —)■ 7 /i+/x“)/rg re¬ 
sult to 8.19 X 10“^ once the empirical central value of 
Pg = r(r 7 — >■ 77) is imposed. The estimates we have 
made agree with the experimental central values at the 


^ We thank P. Masjuan for graciously providing it to us [140]. 

® We thank the authors of Ref. [142] for a high precision version 
of the function P’[{Q‘^) that appears in their Table 5, making it 
suitable for numerical work [140, 142], 
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« 5% level, a disagreement larger than one would naively 
expect from radiative corrections but which, rather, ap¬ 
pears to stem from differences in recent experimental re¬ 
sults for T{r] —)■ e+e“ 7 ) [143, 145] and the world av¬ 
erage [139]. It is worth noting that the ij decay par¬ 
tial widths are measured relative to Fq, and although a 
recent, precision measurement of Fq exists [146], differ¬ 
ences in measurements of Fq from the e“'"e“ collisions 
and electron-nucleus scattering, through use of the Pri- 
makoff effect, impact earlier determinations of the partial 
widths [139]. Irrespective of this, we conclude that our 
framework certainly works sufficiently well for our cur¬ 
rent purpose. Turning to the possibility of an A' inter¬ 
mediate state, we update Eq. (14) to 


F(r7^ 7^' ^ 7£+r) = F2 


X 


2(1 — xYx 


X - XA'f' + 


/3 




\fr^{l,0,x)\^ 



( 20 ) 


where M^o —)• M^, n = 2 mf/M^, and /? = ^X — r^/x. 
Here, too, we can use the narrow width approximation 
to find 


F(?7 ^ 7H' 7£+£ ) = \fr,{l, 0 ,XA')fT{r] 7H') 

xBiA'-^ £+£-), ( 21 ) 

where F (?7 —)■ 7 H') and F(yl' —)■ follow from suit¬ 

able substitutions in Eqs. (17) and (18). 


C. The decay rj —>■ ')A' —>■ 77r"'"7r 



FIG. 3. The branching ratios for B{A' —>■ e''"e“) (solid), 
B{A' —>■ (dotted), and B{A' —>■ (dashed) as 

per Eqs. (18), (23), and (27), assuming that decays to the 
hidden sector do not occur, as a function of the Al mass in 
GeV. 


revealing that Eq. (24) is a very good approximation over 
a large mass range — only the factor (1 — 4M^/s)^/^(l -I- 
2M^/s) is extraneous. A simple form of TV especially 
suitable for y/s < 600 MeV is [148] 


F^{s) 


fp^9p 


s-M^ + m{s) 


(27) 


To study this particular channel we need only compute 
F(A' —)■ 7 r+ 7 r“), replacing £ with tt in Eq. (21). Recalling 
the definition of the tt form factor i7n-(s) in e^e“ —> tt+tt” 
at center-of-mass energy ^/s [147], we note the decay am¬ 
plitude is of form 


M{A'(p) 7r+(pi)7r (pa)) = se£%(pi -p2)pF^{p ^), 

( 22 ) 

so that the decay width is 


1 / \ 

F(A' ^ TT+TT”) = —e^amA' ( 1 - 4^=2 J |EV(to^')P ■ 

12 \ J 

. 

This can be compared to the estimate used earlier in the 
literature [37, 52] 

F(A' —> hadrons) = l^£^amA'F{m\,) (24) 

o 


where the cross-section ratio 

, , _ a{e'^e~{s) —?■ hadrons) 
^ ~ a{e+e-{s)—)■ p+p-) 


(25) 


below TTw threshold evaluates to 

^ 1 (s-4M,^)3A|f,(g)|2 

^ 4(s-4M2)i/2(s + 2M2) ’ 


( 26 ) 


where n(s) = {M^/y£s){p{s)/p{Mp)YTp includes a run¬ 
ning p width, with p{s) = yjsjA — and a nor¬ 

malization as per Ref. [149]. An explicit test of this 
form is shown in Eig. 5 of Ref. [149], for which fp^ = 
0.122(1) GeV^ [150], pp = 5.8, Mp = 769.3MeV and 
Fp = 150 MeV are used. The p is a broad resonance, 
and if its mass and width are determined as per the real- 
axis prescription of Eq. (27), with a s-dependent width, 
as distinct from the location of its pole in the complex 
plane, then the resulting resonance parameters depend 
on the precise form of Ft^{s) [147, 150]. Such a real- 
axis prescription is employed in the determinations aver¬ 
aged by Ref. [139], for which Mp = 775.26(25) MeV and 
Fp = 149.1(8) MeV are reported. We employ this latter 
set of parameters in our analysis here, and the incurred 
differences are negligible. We display the resulting A' 
branching ratios in Fig. 3, and the results are similar, 
though our F(A' —)• tt+tt”) is a bit larger, than those 
determined using the prescription of Eq. (24) [52]. 


1. On r(A' —>■ hadrons) at higher energies 

The proton initial-state radiation mechanism for dark 
gauge boson production can probe A! masses much in 







excess of 600 MeV, so that we now develop an expres¬ 
sion for r(yl' — > 7r+7r“), as well as for r(^' —)■ hadrons), 
in this mass region. To go beyond the prescription of 
Eq. (24) in estimating r(yl' —)■ hadrons) we use Eq. (26) 
at s = to rewrite Eq. (23) as 


r(A' 


-I- —\ ^2 

r^TT ) = -£ arriA' 


1-4- 


Mt 


1/2 




1 + 2-f U(m^0, 

lit' A/ I 


(28) 


which is some 10% smaller than Eq. (24) if niA' = 
300 MeV. For niA' in excess of ttoj threshold, or « 
920MeV, Eq. (28) is simply that of r(yl' —> hadrons), 
implying that B{A' —)• e~^e~) +jB(A' —> +B(A' —)• 

7r+7r“) < 1 as well. 


£. On r(T' —>• invisibles) 

It is entirely possible that a produced dark gauge bo¬ 
son also decays invisibly to the particles of the hidden 
sector, though this possibility has been neglected in early 
studies [37]. Nevertheless, it can readily be incorporated, 
both in this case and throughout, by making the total vis¬ 
ible decay width some fixed fraction of the total width. 
Varying the fraction of decays to the invisible sector con¬ 
sequently changes the exclusion plot. 


D. The decay rj —>■ yp' —>■ yw+Tr 


Now we turn to the possibility of a dark gauge bo¬ 
son, the p', which may or may not be a composite ob¬ 
ject. To anchor our analysis we first compute the rate for 
vip) y(A:)7r+(p2)7r“(pi) in the SM, where we assume, 
in analogy to our earlier work, that the decay vertex is 
mediated by the axial anomaly in the low-energy, chiral 
theory of hadrons [151, 152], where we refer to Ref. [138] 
for an analysis of this assertion within chiral perturba¬ 
tion theory. In particular, we assume that the decay is 
dominated by ry —)■ (7r7r)i=i7 channel with the p playing 
a prominent role. Since p —>■ 7r"''7r“ is a p-wave decay, we 
define the coupling pp as per 

(7r+(p2)7r"(pi)lp°(pp,£p)) = -ppSp ■ ip2-pi) (29) 


to yield the decay amplitude 




* 


{k)sy{pp)kPp: 


fpy9p£p ■ {P2 - Pi) 
s-M^+iU{s) ’ 
(30) 


where s = p^p and n(s) is defined after Eq. (27). Summing 
over the p spin, this becomes 


M = -A^ep,ypae^*{k)pik'’p2 


fp~f9p 


s — + zn(s) 


(31) 


We note that this is just a special case of A4 = 
£pvpa£^*Pik'^P 2 P{s)FT^{s), where P{s) is a polynomial 
of form P{s) = Ap{l -|- as) [153]. The corrections to this 
framework in p —>■ 7r“'"7r“y decay are very small [153-155], 
supporting its use in the analysis here. Evaluating the 
integral over the three-body phase space (with the inte¬ 
gral over solid angle computed in the rest frame of the 
7r+7r“ pair) yields 


r(?y —> y7r''’7r ) 


( 32 ) 


where r = 2M,r/-Mp, x = s/M^, and we have included 
the polynomial of Ref. [153] as Pp{s) = 1 + as with 
a = (1.96(27)stat(2)sys) GeV“^ [153]. Although we could 
determine the constant A^ by simply appealing to the 
empirical rj —>■ y7r“'"7r“ partial width, it is worth noting 
that the axial anomaly can fix this constant as well. That 
is, the low-energy theorem relating y* —)■ to 

7 r° —> 77 [156-158] generalizes in the chiral SU(3)/ limit 
at zero momenta to yield A^ = e/(4-\/37r^/^) [159, 160]. 
Using = 92.2 MeV, this yields = 5.65 GeV“^, 
whereas including SU(3)/ breaking and p — p' mix¬ 
ing gives an additional multiplicative factor of ^ = 
cos^/t/Zs — sin6l-\/2/7r//o [159, 160], which evaluates 
to Z = 1-057 using /g//^ = 1.3, /o//^ = 1.04, and 
9 = —20° [159, 160]. Doing the integral in Eq. (32) using 
Eq. (27), A^ = 5.65GeV■^ and Ttot = 1.31(5) keV [139] 
yields a branching ratio of 6.42%, whereas the experimen¬ 
tal result is B{p —> Tr+Tr^y) = (4.22(8))% [139], a frac¬ 
tional difference of some 30%. Including ^ and an overall 
correction factor to Ap of (1 -l-d) with S = —0.22(4) from 
the assessment of Ref. [153], yields B{p —>■ 7r+7r“y) = 
4.37%, in reasonable agreement with experiment given 
the various uncertainties. To adapt this framework to the 
computation of r(r 7 —)■ yp' —>■ y7r+7r“) we need only re¬ 
place (s-M2)2+n2(s) in Eq. (27) by (s-)2 + m2 T^, 

and multiply by £^. Although B{p' —)■ tt+tt”) could 
plausibly be less than one, we shall assume that it is 
still reasonable to use the narrow width approximation 
of Eq. (15), yielding 


2 ^p' 


r(r 7 ^ yp'^ yTT+TT ) = ^-^M]^\Ap{l + 5)^ 




1 - 


4mA 


3/2 


(33) 


Supposing a p' coupling of the form given in Eq. (29), 
but with strength pp£, we have 


( \ 3/2 

1 ^ j ■ 
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Thus if Tp/ = r(/5' —>■ TT+TT ), we find 


r{r] jp') 





(35) 


In what follows we include the possibility of invisible de¬ 
cays by studying the impact of making B{p' tt+tt”) 
some definite fraction less than one on the parameter ex¬ 
clusion plot. 


pairs and are subsequently detected in the pair spectrom¬ 
eter. 

We have been conservative in the following simulations. 
We have only considered events where the background is 
known to be low, namely after all or most of the shield. 
For the special case of muon pair detection, we have used 
part of the last meter of the Fe shield as part of the 
fiducial region since the Fe absorber has little effect on 
the muons. As more is known about the experiment and 
if more trigger optimization can be implemented, it may 
be possible to extend the fiducial region further upstream 
in the beam dump. However, such studies are beyond the 
scope of this work. 


IV. EXPERIMENTAL PROSPECTS 

Already, there are a large number of searches for dark 
photons and a number of planned searches [37, 40-52, 54- 
65]. 

We illustrate here a possible new search for dark pho¬ 
tons at the FNAL E906 (SeaQuest) experiment at Fer- 
milab. This experiment was designed to perform a mea¬ 
surement of Drell-Yan processes for a proton beam inci¬ 
dent on stationary targets. The experiment is comprised 
of an intense (5 x 10^^ protons/spill at 1 spill per min) 
120 GeV proton beam incident on a target. Less than 
a meter downstream of the target is a 5-m long solid 
iron magnet that not only serves to begin the analysis of 
the muon pair from the Drell-Yan process, but also as a 
beam dump. This beam dump tremendously attenuates 
all hadrons and most charged particles from the target 
with the exception of muons. Following the solid iron 
magnet, there are two stations of scintillator hodoscopes 
and drift chambers that can define the vertex of either 
a prompt Drell-Yan event in the target or beam dump 
or from the decay of a prompt or displaced vertex from 
a dark photon. The first detector station is followed by 
an air gap magnet and two additional detector stations. 
This system serves as the pair spectrometer and can mea¬ 
sure the mass of an event with a relative mass resolution 
that varies between 2.5-6% depending on where the event 
occurs in the solid iron magnet. A 1-m thick block of iron 
was placed just upstream of the fourth detector station 
and serves to filter most charged particles except muons. 
The fourth detector substation comprises a layer of pro¬ 
portional tubes and scintillator paddles that serve as a 
muon identification system. A schematic overview of the 
SeaQuest spectrometer is shown in Fig. 1 [36]. 

Clearly, experiment E906 has all the basic elements of 
a “shining-through-the-wall” beam stop experiment to 
search for dark photons. The basic elements consist of a 
high-energy proton beam incident on a fixed target and 
a pair spectrometer which detects lepton pairs from the 
decay of the dark photon, where the spectrometer is well- 
shielded from the target in order to minimize background 
events from ordinary SM processes. Consequently, the 
dark photons, being weakly interacting, can travel un¬ 
scathed through the shield where they decay into lepton 


A. rj meson decay to charged lepton pairs 

Here we illustrate the regions of sensitivity in parame¬ 
ter space of e and dark photon mass, toa', for SeaQuest 
for four different processes that could be used to search 
for dark photons. Fig. 2 shows the projected regions 
of sensitivity for SeaQuest along with areas excluded by 
electron and proton beam dump experiments and pro¬ 
posed experiments. The simulation assumes that we have 
a 200-day experiment with an overall efficiency of about 
30% which approximates the conditions experienced in 
the experiment thus far. We investigated dark photon 
production from radiative rj and decays and proton 
bremsstrahlung. The rj and 7r° yields per proton have 
been estimated from GEANT4 simulations that record 
the energy and transverse momentum spectra of rj and 
7r° mesons produced (shown in Fig. 4 and Fig. 5) when 
6.2 million 120 GeV protons are incident on a 4.75-m 
long Fe beam dump. These energy and px distributions 
were used in the calculation of the estimation of the dark 
photon yield. 



FIG. 4. The rj (solid) and 7r° (dotted) yield/proton as a func¬ 
tion of the energy of the particles as obtained from GEANT4 
Monte Garlo simulations of 120 GeV protons interacting with 
the Fe beam dump. 
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FIG. 5. The 77 (solid) and 7 r° (dotted) yield/proton as a func¬ 
tion of the transverse momentum pr obtained from GEANT4 
Monte Garlo simulations of 120 GeV protons interacting with 
the Fe beam dump. 



FIG. 6 . Preliminary projections for an 77 meson decaying to 
a photon (7) and a dark photon, denoted by A'. The A' can 
subsequently decay into a p,~ pair (solid), or an e^e“ pair 
(dashed). Here decays to the hidden sector do not occur. The 
dotted curve represents the projected exclusion limit for tt® 
decay. 


In this projection, we assume that the fiducial region 
for detecting the decay of the dark photons to a muon 
pair begins 4-m downstream of the front edge of the beam 
stop and extends 0.85 m downstream. This hducial re¬ 
gion ensures that all decay muons must travel through at 
least 0.15 m of FMAG, which would impart an effective 
Pt kick of about 0.08 GeV. This ensures that the muon 
pair would be well separated before reaching the first de¬ 
tector station. For the decay of a e“'"e“ pair or in later 
discussion, a pion pair, the fiducial region is only 0.95 m 
and entirely outside the Fe magnet. For the case of e+e” 



FIG. 7. Here the A' is produced by 77 decay. The A' sub¬ 
sequently decays to a muon pair. The preliminary projected 
exclusion regions for SeaQuest for three branching ratios to 
the dark sector between 0 and 90% are shown. Specifically, 
B{A' —>■ invisibles = 0,0.70,0.90, which are shown as solid, 
dotted, and dashed curves, respectively. 


or 7r“'"7r“ pairs, it might be necessary to add a small air 
gap magnet between the downstream end of FMAG and 
the first detector station to ensure a small pp kick. For 
the purpose of simulation, it was assumed that the nec¬ 
essary separation of the pairs was met. To be specific, 
the expression used for the A' decays to lepton pairs in 
the fiducial length, Zfid, after the A! traverses the length 
of the beam dump, ^dump, is given by 


iVdec = iVoS(A' ^ £+r)exp 

\ cta' \pa'\J 


1 — exp ( — 


^fid ^A' 

cta' \pa’ 


(36) 


where pa' is the A! momentum in the laboratory frame, 
and ta' is the A lifetime in the A rest frame. 

Two different approaches were used to determine the 
acceptance of the spectrometer for a dark photon. They 
both assumed a dark photon of various masses decay¬ 
ing at different vertex positions in the beam stop. The 
first method makes use of the full spectrometer geom¬ 
etry in a GEANT4 (version 4.9.6.p03 [161]) simulation. 
The second method uses the spectrometer geometry in a 
simplihed situation where multiple scattering and energy 
loss of the decay particles are put in explicitly. The kine¬ 
matic variables and the dimuon yields were compared, 
and there is good agreement between the two approaches. 
Gurrently, our acceptance is limited by the decay length 
of the dark photon. Once track reconstruction is opti¬ 
mized and a detailed simulation of the background events 
are achieved, it may be possible to probe larger values in 
£ parameter space. 

The projected regions of sensitivity based on this simu- 
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Mp. (MeV) 


FIG. 8. Preliminary projections for an r; meson decaying 
to a photon and a dark gauge boson p' that couples to 
quarks, with a probability governed by the square of the cou¬ 
pling constant between light and dark sectors of with the 
probability governed by Eq. (35). The p' can subsequently 
decay into a charged pion pair with a decay lifetime pro¬ 
portional to governed by Eq. (34). The p' is permit¬ 
ted to decay to the dark sector with a branching ratio of 
B{p' — >■ invisibles) = 0,0.80,0.95, which are shown as solid, 
dotted, and dashed curves, respectively. 


lation are indicated in Fig. 6 . The contours shown in the 
figure indicate 10 -event contours, which should be ade¬ 
quate to set 95% conhdence levels on exclusion of dark 
photons. For example, if we assume a signal of 10 events 
and a detected number of events of 10 with a background 
of 7 events, the confidence level [162] that the true num¬ 
ber of events is less than 10 is 94.6%. Three cases are 
displayed: ( 1 ) detection of a muon pair from 77 decay to 
a dark photon and subsequent A' decay to a muon pair, 
( 2 ) detection of an e+e“ pair from rj decay and subse¬ 
quent A' decay, and (3) detection of e+e“ decay from 
neutral pion decay and subsequent A' decay. It is as¬ 
sumed that there are no decays to the invisible sector for 
the projections of Fig. 6 . 

Many projections and exclusion plots assume that the 
decay to the visible sector is 100%. In Fig. 7 we explore 
the case where A' decay to the invisible sector is also 
allowed to determine its impact on the ability to limit 
the A'. The figure shows the specific case of 77 decay 
to 7 %', followed by A' —■ Furthermore, it is 
assumed that the branching ratio to the hidden sector, 
B{A' invisibles, varies from 0 to 90%. Clearly, if decay 
to the hidden sector is permitted, the region excluded by 
the experiments dramatically shrinks. 


B. rj meson decay to a charged pion pair 

With relatively small modifications, it may be possible 
for the pair spectrometer associated with the SeaQuest 
experiment to detect a charged pion pair. Although the 
detectors in this spectrometer presently do not have the 
capability to distinguish between an e“'"e“ pair and a pion 
pair, it may not be necessary in a first stage of the ex¬ 
periment. If a mass peak is observed from a displaced 
vertex in the SeaQuest experiment, then the necessary 
particle identification could be added to the detector to 
determine whether the pair is leptonic or hadronic. If 
we consider 77 decay to a photon and a dark p', and the 
dark p' subsequently decays to a pion pair, the expected 
sensitivity plot is given in Fig. 8 , again showing 10-event 
contours based on Eqs. (35) and (34). The effect of the 
p' decay to the invisible sector has also been explored in 
the figure. 

V. SUMMARY AND OUTLOOK 

In this paper we have explored the discovery prospects 
of the SeaQuest experiment, E906 at Fermilab, operated 
in proton beam dump mode. We have explored the abil¬ 
ity to which it can explore new regions of dark photon 
parameter space, and how, if pion pair detection were 
made possible, it can search for an entirely new kind of 
hidden sector, one probed through the hidden strong in¬ 
teractions of quarks. The latter scenario has never been 
explored in beam dump experiments, and the possibil¬ 
ity of strongly coupled hidden sectors at sub-GeV mass 
scales is poorly constrained by other considerations [ 110 ]. 
We emphasize that the SeaQuest projectionss in this pa¬ 
per have been computed assuming one year of running 
with a non-optimized detector. Increasing the exposure 
could lower the epsilon range considerably, and, in the 
case of the proton bremsstrahlung mechanism, permit 
sensitivity to higher-mass gauge bosons as well. 

Moreover, the prospect of polarized proton beams [126] 
opens a unique window of sensitivity to light Z' gauge 
bosons, to the Zd scenario of Refs. [67, 84, 93, 94] and, in 
the event of tt+tt” detection, to the leptophobic, light Z' 
proposed in Ref. [110]. Increasing either the proton beam 
energy or the length of the fiducial region in the iron 
beam stop can also extend the sensitivity of SeaQuest to 
heretofore unexplored regions of hidden-sector parameter 
space. 
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